New constraints on parametrised modified gravity 
from correlations of the CMB with large scale structure 



o : 

(N . 
^ ■ 

<: 

in ■ 



u: 
6 ■ 

(N ■ 

> ■ 

in : 
o 

0^ ■ 

o : 

o\ . 

o . 



Tommaso Giannantonio,^ Matteo Martinelli,^ Alessandra Silvestri,^ and Alessandro Melchiorri^ 

^ Argelander-Institut fiir Astronomie der Universitdt Bonn, Auf dem Hiigel 11, D-53121 Bonn, Germany 
^Universitd degli studi di Roma "La Sapienza", P.le Aide Mora 5, 1-00185 Roma, Italy 
'^Kavli Institute for Astrophysics and Space Research, MIT, Cambridge, MA 02139 USA 

We study the effects of modified theories of gravity on the cosmic microwave background (CMB) 
anisotropics power spectrum, and in particular on its large scales, where the integrated Sachs- 
Wolfe (ISW) effect is important. Starting with a general parametrisation, we then specialise to 
f[R) theories and theories with Yukawa-type interactions between dark matter particles. In these 
models, the evolution of the metric potentials is altered, and the contribution to the ISW effect 
can differ significantly from that in the standard model of cosmology. We proceed to compare these 
predictions with observational data for the CMB and the ISW, performing a full Monte Carlo Markov 
chain (MCMC) analysis. In the case of f{R) theories, the result is an upper limit on the lengthscale 
associated to the extra scalar degree of freedom characterising these theories. With the addition of 
data from the Hubble diagram of Type la supernovae, we obtain an upper limit on the lengthscale of 
the theory of Bq < 0.4, or correspondingly Ai < 1900Mpc//i at 95% c.l. improving previous CMB 
constraints. For Yukawa-type models we get a bound on the coupling 0.75 < /3i < 1.25 at the 95% 
c.l. We also discuss the implications of the assumed priors on the estimation of modified gravity 
parameters, showing that a marginally less conservative choice improves the f{R) constraints to 
Ai < 1400Mpc//i, corresponding to Bo < 0.2 at 95% c.l. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 



I. INTRODUCTION 

The observed acceleration of the cosmic expansion 
poses a challenge for modern cosmology. Standard gen- 
eral relativity (GR), applied to an expanding universe 
filled with radiation and cold dark matter, does not fit 
the data unless one invokes an additional component, 
either a cosmological constant, A, or a dynamical dark 
energy field [l| . The former case corresponds to the well 
known standard model of cosmology, or ACDM. An al- 
ternative approach to the phenomenon of cosmic accel- 
eration consists of modifying the laws of gravity on large 
scales, in order to allow for self-accelerating solutions. 
Well-known examples of such theories are the f{R) mod- 
els or more general scalar-tensor theories [7|-[l0|. 
the Dvali-Gabadadze-Porrati (DGP) model [Til. [l2| . and 
its further extensions such as Degravitation [l3|. 

Such theories need to fit the observed data as well as 
or better than the ACDM baseline model. This require- 
ment is twofold: at the background level, the predicted 
expansion history has to be coherent with distance mea- 
surements such as those from Supernovae (SNe), the po- 
sition of the CMB peaks, and baryon acoustic oscillations 
(BAG); on the other hand, the modified predictions for 
structure formation need to pass the test of comparison 
with the observed large scale structure (LSS) of the Uni- 
verse. It has been shown for several cases — for instance 
for f{R) [14- 18] and DGP [H-Ill — that the structure 
formation test can dramatically reduce the parameter 
space left unconstrained from the background test, and 
it is therefore instrumental to the quest for a truly viable 
model. 

Structure formation in the Universe is driven by grav- 



ity through its potentials, which at the perturbative level 
in Newtonian gauge are described by two fields *I'(a:^) and 
^{Xfj), defined as the perturbations to the time-time and 
space-space parts of the metric respectively. A powerful 
test of modified gravity is offered by measuring the past 
evolution of these potentials to compare it with the pre- 
diction of a given theory. In the standard ACDM model 
the potentials are expected to be identical, to remain 
constant during matter domination and to decay at late 
times, as a consequence of the ongoing transition to the 
dark energy phase. This is not necessarily the case for 
other theories. Metric perturbations are not viable as- 
tronomical observables, but we may reconstruct them 
from observational data which are directly dependent on 
them [2^, [2^ . The most useful data for this purpose are 
galaxy and cluster counts, weak lensing and the CMB. 

The CMB is an almost perfectly isotropic black body 
radiation that has been generated at the epoch of hy- 
drogen recombination and has undergone free streaming 
since then. Nonetheless, this free streaming may be al- 
tered if the CMB photons encounter potential wells which 
evolve in time. In this case, a non-zero net energy is 
gained (or lost) by the photons if the potential is be- 
coming shallower (or deeper): this phenomenon, known 
as the integrated Sachs- Wolfe (ISW) effect [131 , leaves 
an imprint on the larger scales of the CMB spectrum. 
In the standard ACDM theory, this effect is expected to 
be generated at late times as a consequence of the po- 
tential decay when the background starts accelerating. 
In models of modified gravity, such as f{R) theories, the 
magnitude of this late effect can be altered, as well as po- 
tentials can evolve during matter domination, therefore 
generating an ISW effect at earlier redshifts [l3, [13 ■ 
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The ISW effect can be measured by cross-correlating 
the CMB with tracers of the large scale structure 
(LSS) and it has been detected using several differ- 
ent data sets and was used to constrain cosmology [29'- 
lilj . The strongest (4.5 a) detection to date [42] has been 
obtained by combining multiple data sets at different red- 
shifts, thus exploring the evolution of the potentials in 
time, and is fully consistent with the ACDM picture. 

Here we focus on f{R) theories and models with a 
Yukawa-type dark matter interaction, and use the ISW 
data by [42|, in conjunction with the CMB, to test the 
general parametrisation of modified gravity by exploring 
the parameter space with a Monte Carlo Markov chain 
(MCMC) technique. 

The plan of the paper is as follows. In Section HIl first 
we introduce the parametrisation used to describe depar- 
tures from GR, and then specialise to the case of f{R) 
and Yukawa-type theories. Then we review the method 
of the analysis and the data used in Section IIIII In Sec- 
tion |IV] we give details of our results and in Section IVl we 
discuss the effects of the assumed priors on the parameter 
estimation, before concluding in Section IVTl 



II. PARAMETRISED MODIFIED GRAVITY 

Here we describe the formalism we use to parametrise 
departures from general relativity. 



A. Background expansion 

In our analysis we fix the background to that of the 
ACDM model of cosmology. The reasons for this choice 
are multiple; ACDM is currently the best fit to avail- 
able data and popular models of modified gravity, e.g. 
f{R), closely mimic ACDM at the background level with 
differences which are typically smaller than the precision 
achievable with geometric tests [isj . The most significant 
departures happen at the level of growth of structure and, 
by restricting ourselves to ACDM backgrounds, we can 
isolate them. 



B. Structure formation 

In models of modified gravity we expect departures 
from the standard growth of structure, even when the 
expansion history matches exactly the ACDM one. In 
general, the rate of clustering of dark matter, as well as 
the evolution of the metric potentials, is changed and can 
be scale-dependent. Moreover, typically there might be 
an effective anisotropic stress introduced by the modifica- 
tions and the two potentials, $ and are not necessarily 
equal as it was the case for ACDM [H, [H, |T3, M- 
[22l . |44| . Here we focus on the effect of the modi- 
fied evolution of the potential, $ -t- ^t, on the CMB. 



In order to evolve the potentials we employ the MG- 
CAMB code developed in p.8] (and publicly available 
at |http : //www . sf u . ca/~gza5/MGCAMB . html|) to evalu- 
ate the growth of perturbations in models of modified 
gravity. In this code, the energy-momentum equations 
remain the standard ones, and the modifications to the 
Poisson and anisotropy equations are encoded in two 
functions /Lt(a, k) and 7(0, k) defined by 



2Af|, 



fi{a,k)pA , 



— =7(a,fc). 



(1) 
(2) 



where pA = pS + 3^(p -f- P)v is the comoving density 
perturbation. 

We will consider theories in which the modifications in- 
troduce an effective scalar degree of freedom (d.o.f.), and 
therefore a characteristic lengthscale. A typical action 
for such theories would be 



S 



M2 



-a,(4>}/Mf 



1 



g^" (V^(/.) (V,0) - y(0) 



(3) 



where (j) represents the scalar d.o.f., Xi is the r mat- 
ter field and Q;,i(0) is the coupling of to Xi- We will 
limit ourselves to cases in which the coupling is a linear 
function of the scalar field, i.e. ai{(j)) oc 0. 

For the theories described by the action in Eq. ([3]) , the 
functions p and 7 can be well represented by the following 
parametrisation introduced by (45l | (and used in [Tsl )) 



p{a,k) 



l{a,k) 



1 + PiXl Pa' 
1 + Xj k^a'' 

1 + A2 k^a' 



(4) 
(5) 



where the parameters /3i can be thought of as dimcn- 
sionless couplings, as dimensionful lengthscales and s 
is determined by the time evolution of the characteristic 
lengthscale of the theory, i.e. the mass of the scalar d.o.f. 
As shown in [l^ , in the case of scalar-tensor theories the 
parameters A^} are related in the following way 



/3i 



A? 



'Af 



(6) 



and 1 < s < 4. 



1. f{R) theories 

These theories are a subclass of the models described 
by the action of Eq. ([3]), corresponding to the case of a 
universal fixed coupling = (f>, i.e. /3i — 4/3 [il]. 

Moreover, f{R) models that closely mimic ACDM cor- 
respond to s ~ 4 fl8|. Therefore the number of free 
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parameters in Eqs. ^ and ([5]) can be reduced to one, 
e.g. the lengthscale Ai. 

The parametrisation in Eq. Q effectively neglects a 
factor representing the rescaling of the Newton's con- 
stant already at the background level due to the mod- 
ifications (e.g. (1 -I- in f{R) theories). Such a 
factor is very close to unity in models that satisfy local 
tests of gravity [i^ and, as such, it can be neglected. 
However, when studying the f{R) case, we need to in- 
clude it to get a more precise estimate of the ISW effect 
for the MCMC analysis; therefore we use the following 
extension of Eq. (|4]) 

1 1 + -Xj fc^a^ 

" l-1.4.10-8|Ai|2a3 T+Affc^ ' 

where the difference with Eq. (|4]) is the multiplica- 
tive factor parametrising (1 + fR)~^ which can be ex- 
pressed in terms of Ai [4J, |43|. Even with this extended 
parametrisation, we are left with a single free param- 
eter, the lengthscale Ai. This parameter can be eas- 
ily related to the value of the mass scale of the scalar 
d.o.f. introduced by the addition of the f{R) term to 
the Einstein-Hilbert action. This scalar d.o.f. is rep- 
resented by the function = df/dR, also dubbed the 
scalaron, and Ai corresponds to its mass scale today, i.e. 
Ai = In [11 Ea] this family of f{R) models was 

labelled by the parameter Bq which is related to Ai via 
A^ = Bq c? /{2Hq). We will present our results mainly in 
terms of Bq to facilitate the comparison. 



2. Yukawa-type dark matter interaction 

Scalar-tensor theories are an example of models with 
interaction between dark energy and dark matter. The 
parametrisation in Eqs. ([7]) and ([5]) is well suited also for 
models where dark matter (DM) particles self-interact a 
la Yukawa. Such models can be described by the presence 
of an additional mediator, represented by a scalar d.o.f., 
like the one in the action of Eq. ([3]) , with a time-evolving 
mass and a coupling to DM particles which is now a free 
parameter of the model. In this case the only coupling 
of interest in Eq. ([3]) would be the one to DM, adm 7^ 0. 
The coupling of this scalar d.o.f. to other matter fields, 
such as baryons, does not need to be the same of the one 
to DM; in particular, given the stringent constraints on 
such a coupling, it can be considered negligible. 

In our parametrisation, the additional mediating 
force is characterised by the coupling /3i (related to 
ttdm), the lengthscale Ai and its time evolution s. The 
remaining parameters, (32 and A2, can be related to /?i 
and Ai as in the scalar-tensor case of Eq. ([5]). In this 
case there is no need of any additional factor in the 
expression for /i, since it takes already into full account 
the effect of a Yukawa coupling on the clustering of 
dark matter, (in other words the Yukawa interaction 
alters the potential produced by a dark matter particle. 



but does not alter the background Newton's constant). 
Therefore, when analyzing models with Yukawa-type 
dark matter interaction we will use Eqs. ([7]) and ([5]) with 
free parameters {/3i, Ai, s}. 

Armed with the expressions for and 7 to feed into 
MGCAMB, we can proceed to evaluate our observables. 
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FIG. 1: Top: theoretical predictions for a family of f{R) 
theories compared with our ISW data [i^l measuring the an- 
gular CCF between the CMB and six galaxy catalogues. The 
model with Bo = is equivalent to ACDM , while increas- 
ing departures from GR produce negative cross-correlations. 
Bottom: the same for a family of Yukawa-like theories with 
fixed Bo ~ 2. In this case non-unitary coupling generate a 
redshift evolution of the signal. 
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Parameter 


Explanation 


range (min, max) 




Primary parameters 


f{R) 1 Yukawa-type 


r 

Asz 

Us 


physical baryon density; ujb ~ /I'^f^b 

physical CDM density; luc = h^Qc 

sound horizon angle; i9» = 100 ■ rs{z,) /Da{z,) 

optical depth to reionisation 

curvature density; Q,k = 1 — fitot 

As is the scalar primordial amplitude at fepivot = 0.05 Mpc~^ 
amplitude of the SZ template for WMAP and ACBAR 
spectral index of primordial perturbations; Us — 1 = d\nP/d\nk 


(0.005,0.100) 
(0.01,0.99) 
(0.5, 10.0) 
(0.01,0.80) 


(2.7,4.0) 

(0,2) 
(0.5, 1.5) 


Bo 

s 


present lengthscale of the theory (in units of the horizon scale) 
coupling 

time evolution of the scalaron mass 


(0,6) 
4/3 
4 


(0,6) 
(0.001,2) 
(1,4) 




Derived parameters 






Ho 
h 

Ha 


Hubble parameter [km/s/Mpc]; calculated from uji,, ujc, i9, and fix 
h = Ho/(WO km/s/Mpc) 

matter density parameter; Qrn ~ {i^b + (^c)/h'^ 
vacuum energy density parameter; Qa = 1 — — 


tophat (40, 100) 
(0.40, 1.00) 


Ai 

/52 


lengthscale of the theory; Ai = Boc^/{2Ho) 
second lengthscale of the theory; A2 — /SiAj 
anisotropy parameter; 2//3i — 1 


(0,5- 100 
(0,7- 10^ 
1/2 


(0,5 • 10') 
(0,10«) 
(0, 2000) 



TABLE I: Our primary MCMC sampling parameters and some important derived parameters. The modified gravity parameters 
are separated by a horizontal line. 



III. ANALYSIS 
A. Method 

We run Monte Carlo Markov chains using the stan- 
dard Cosmomc package [ioj . which is based on the 
CAMB CMB code [50| and is publicly available at 
[http : //www, cosmolog ist . info/ cosinoinc[ As it is well 
known, this technique consists of using a Metropolis- 
Hastings algorithm to efficiently sample the multi- 
dimensional parameter space in which our model lives. 
The algorithm is based on multiple chains, which are 
started from different random initial points, and then 
evolved based on a prior probability distribution of the 
parameters, assumed flat in the chosen range. Each new 
step will be accepted if its likelihood improves the pre- 
vious likelihood, weighted with some probability. The 
final posterior probability distribution is then directly 
obtained from the density distribution in the parameter 
space. 

To study our modified gravity (MG) model we run five 
MCMC chains. We use the parameters described in Table 
|T) this is based onto the usual set of ACDM parameters 
(note that the angle to the last scattering surface z?* is 
used instead of the Hubble parameter Hq since this re- 
duces the degeneracies with the other parameters). To 
account for f{R) theories we have to add one extra pa- 
rameter: here we chose to use to use Bq as a primary 
parameter. The other MG parameters are then fixed: 
\l = c-'Bo/{2H^), Ai = /3i = 4/3, /32 = 1/2 and 

s = 4. In theories with Yukawa-type dark matter inter- 
action we have two additional free parameters, i.e. the 
coupling /?! and s as described above. 

We shall discuss the effects of the parameter choice on 



the priors and on the results in more details in Section 

El 

B. Data 

1. CMB 

For the CMB we use the publicly released WMAP 5 
years data [5lj for the temperature and polarisation (TT, 
TE and EE) power spectra of the perturbations. In ad- 
dition, we used flat top hat priors on the age of the Uni- 
verse, to E [10, 20]Gyrs and on the Hubble parameter to- 
day Hq e [40, 100]km/s/Mpc. We include CMB lensing 
in the analysis. 

2. ISW 

To break the degeneracy which remains at the back- 
ground level between GR and MG, we use the ISW data 
by [42]. These data were obtained by cross-correlating 
the WMAP maps of the CMB with six galaxy data sets 
in different bands (2MASS, SDSS main galaxies, LRGs 
and QSOs, NVSS, HEAO). The data span different red- 
shift ranges from z = 0.1toz = 1.5, thus allowing us to 
study the evolution of gravity, through the history of the 
decay of its potentials, in a tomographic way. 

There are 78 actual data points, consisting of the an- 
gular cross-correlation functions (CCFs) in real space, 
binned at 12 angles between and 12 deg for each of 
the catalogues. The full covariance matrix is highly non- 
diagonal both between the different angular bins and be- 
tween the catalogues, due to overlaps in redshift and sky 
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coverage. For each MC model, the hkehhood contribu- 
tion due to the ISW data is calculated as follows: first the 
theoretical matter and matter-temperature power spec- 
tra are calculated through a full Boltzmann in- 
tegration inside CAMB, then a Legendre transformation 
yields the CCFs and the matter auto-correlation func- 
tions (ACFs) at the relevant scales. Finally, the bias 
parameters, assumed constant for each catalogue, are re- 
calculated for each model by forcing the ACFs to match 
the observations. 

We show in Fig. [T] the CCF data points, overlapped 
with the theoretical predictions for a family of f{R) and 
Yukawa-like theories. We can see that in the first case 
a departure from GR produces an unobserved negative 
signal, while in the second a coupling different from unity 
causes an equally unobserved redshift evolution of the 
effect. 



adding these constraints, the parameter space allowed 
from the background expansion history becomes nar- 
rower, and thus so become the constraints on the mod- 
ification of gravity, mainly due to the degeneracy of the 
theory wavelength Ai with other parameters such as ^m, 
as can be seen in the 2D contour plot of Fig. [3] (top panel). 
We can see the one-dimensional marginalised likelihood 
curves for the models in the left panel of Fig. [5] (black 
solid lines). The resulting constraint from the combina- 
tion of CMB+ISW-f SN is 

Bo < 0.4 or Ai < 1900Mpc//i @95%c.l. (9) 

We have also tested how the results can improve by 
adding the latest additional priors on the value of the 
Hubble constant Hq 53]; however, since the background 
is already well constrained by the intersection of CMB 
and SNe data, the improvement is marginal (at the per- 
cent level in the value of the upper limit). 



To further constrain the background expansion history 
of the Universe, we study the effect of including the con- 
straints from the Hubble diagram of distant Type la Su- 
pernovae (SNe). In particular we use the Union SN com- 
pilation by [53, which consists of 414 SNe drawn from 
13 independent data sets plus 8 newly discovered SNe at 
low redshift, all reanalysed in a consistent way. 



IV. RESULTS & CONSTRAINTS 

A. f{R) theories 

In order to be conservative, we first perform an MCMC 
analysis using the constraints from CMB+ISW only, as it 
is in general good practice to add data sets gradually, to 
control that there are no strong tensions between them, 
which would produce artificially tight posteriors. We 
show in the left panel of Fig. [2] (red dashed lines) our con- 
straints from this run. Here we see the one-dimensional 
marginalised likelihood distributions of the cosniological 
parameters. It can be seen that the usual ACDM pa- 
rameters have likelihoods peaked around their standard 
values, while the extra parameter Bq has an upper limit, 
so that we find 

Bo < 0.5 or Ai < 2000Mpc//i @95%c.l. (8) 

We have checked the convergence of this and the follow- 
ing results from the R — 1 statistic between the chains, 
requiring always R — 1 < 0.02 and over 3000 estimated 
independent samples. This result, which has been made 
possible by a full combined CMB-ISW analysis data, im- 
proves the previous constraints of Bq < 1 by [i^, ob- 
tained only by considering models with a positive ISW- 
matter correlation signal. 

We then perform a new analysis by adding the con- 
straints from the Union supernovae catalogue by (52j . By 



B. Yukawa— type dark matter interaction 

Finally we perform an analysis with three parameters 
to investigate theories in which dark matter particles in- 
teract via a Yukawa-type force. In this case, we directly 
include all data sets, i.e. CMB ISW and SNe, since we 
expect the marginalised constraints to be weaker due to 
the additional parameters. We can see the result of the 
ID posterior likelihood in the right panel of Fig. [5] In this 
case, the upper limit on Bq (or A^) is weakened and be- 
comes uninteresting due to additional degeneracies, while 
we find a constraint on the coupling parameter 



0.75 <l3i< 1.25 @95%c.l. 



(10) 



This limit is enhanced by the lack of strong redshift evo- 
lution of the ISW signal. The result is comparable with 
the previous constraints by [HJ] , although obtained under 
different assumptions. We can not find any significant 
bound on the additional parameter s. We show in the 
bottom panels of Fig. [3] the 2D likelihood contours for 
rim, Bo, f^m,/?i and /3i,Bo. 

From the last panel of Fig. |3]we can better understand 
why the constraint on Bq is now weakened. In the f{R) 
case, the coupling parameter is fixed to /3i = 4/3, for 
which high values of Bq are disfavoured. When the cou- 
pling is set free in the Yukawa- type theories, its preferred 
values are around /3i ~ 1, for which we can see that Bq 
is unconstrained. 



V. ON THE PARAMETRISATIONS AND THE 
BAYESIAN PRIORS 

In Bayesian theory, the posterior probability of a model 
M with N paraniters given the data D is obtained as 



v{e) 



Z{M) 



(11) 
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FIG. 2: Left: posterior likelihood distributions for the model parameters for the f{R) case, using the combined CMB and 
ISW data (red, dashed lines) and adding also the SNe data (black, solid). We can see that the usual ACDM parameters peak 
around the concordance values, while the extra parameter Bo has an upper limit. The SNe tighten the constraints by reducing 
the degeneracy between ilm and Bq. Right: the same for the Yukawa-type case. Here we only show the full CMB+ISW+SN 
result. The upper limit on Bq disappears in this case due to the removal of the corrective factor from Eq. (O and the additional 
degeneracies, and we can not constrain s either, but a value of the coupling /3i close to unity is required to fit the data. 
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FIG. 3: 2D posterior likelihood distributions; the 68% and 
95% probability contours are shown. In the top panel we can 
see the parameters Q,m,Bo for the f{R) case, using the com- 
bined CMB and ISW data (red, dashed lines) and adding also 
the SNe data (black, solid). The SNe tighten the constraint 
by reducing the degeneracy between flm and Bo • In the bot- 
tom panels, we show the Yukawa-type case for the parameters 
nm,Bo (left), (centre), and Pi, Bo (right). 



where £(0) is the likelihood function, 11(0) the prior 
and Z{M) is the Bayesian evidence for the model, which 
we can here consider constant. The marginahsed poste- 
rior of a single parameter ?9 is then obtained integrating 
over the other parameters: 



V{d)(x y"£(0)n(0)d^-i0. 



We normally assume the priors to be flat over the sam- 
pled interval, so that they can be ignored in the marginal- 
isation process. In this case, the posterior is simply 
mapped by the likelihood, upon which assumption the 
MCMC technique is based. However, if we change one 
parameter as ■)? — >■ Cl*^); we have that a prior which was 
flat in -d will not in general be flat in (. In other words, 
changing parametrisation modifies the prior from n('i9) 
to 



n(C) = n(z?)f , 



(13) 



(12) 



which corresponds to a different weighting of the likeli- 
hood in Eq. (jl2p : in higher dimensions, the derivative is 
replaced by the Jacobian of the transformation, as de- 
scribed for example in ^5|. If we now run an MCMC 
chain using the new parameter ^ as a primary param- 
eter, we will in practice force the prior to be flat, 
meaning that the prior on the old parameter n(??) will be 
tilted by a factor dC,/dd. If a prior is tilted, it will cause a 
downweighting of the parameter region in one direction, 
thus affecting the results. All this is well known, and it 
has been discussed by (49j . 

Clearly, if the theory is well constrained by the data, 
then the interval of integration will be small, and the 
prior will be approximately flat in any parametrisation, 
yielding consistent posteriors. However, this does not 
apply to our case: since the current data are not yet 
strongly constraining MG, we have found that the choice 
of priors does indeed have strong effects. In this case, we 
read in [49|i that 

" ...for the results of the parameter estimation 
to be meaningful it is essential that the priors 
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on the base set are well justified..." 

This leaves us with the problem of deciding which of the 
many possible MG parameters (Bo, Ai, A^, LogA^, ...) is 
most meaningful or well justified to be assumed having 
flat priors. The conclusion seems uncertain since any 
modification of gravity is currently based on specula- 
tion without experimental backing, and it is thus hard 
to decide which parametrisation ought to be preferred 
on physical grounds. Since there is no evidence for a 
modification of GR, we have then decided to make the 
most conservative choice by presenting the parametrisa- 
tion which yields the weakest constraints on this modifi- 
cation. 

Nevertheless, at least in the simplest f{R) case with 
one extra parameter, we can forecast the effect of the 
different options by using Eq. (|13p . If we take the pa- 
rameter Ai as our baseline choice, we can predict that by 
switching to fiat priors in the other parameters we will 
have the effects summarised in Table |TTJ any transforma- 
tion to a parameter C, with (PC,/d{\iY > will tilt the 
priors — and bias the posteriors — towards higher values 
of Ai; the opposite for d^(/d(Ai)^ < 0. 

Finally we compare these predictions with the results 
we infer on the baseline parameter Ai by running MCMC 
chains using different primary parameters; we find a good 
agreement with the prediction, and in particular we con- 
firm that our choice of Bq (or equivalently A^ ) is the most 
conservative one. If we decide that the wavelength Ai 
should be considered the most physically justified quan- 
tity, we then obtain the following stricter bounds in the 
f{R) case using CMB+SN+ISW: 

Ai < 1400Mpc//i or Bo < 0.2 @95%c.l. (14) 



formed a full MCMC analysis and found some new con- 
straints on this parameter: in the f{R) case, we have 
obtained Bq < 0.4 or Ai < 1900Mpc//i at the 95% c.l. 
when using combined data from the ISW, the CMB and 
SNe. This puts strong bounds on the parameter space of 
such a theory in light of current data, therefore further 
reducing the possibilities for significant infrared modifi- 
cations of general relativity. 

We have highlighted how the choice of the parametri- 
sation, through its effect on the Bayesian priors, can be 
a very sensitive issue in the estimation of non-standard 
cosmological parameters. We have shown that by using 
a marginally less conservative prior we can reduce the 
posterior constraint to Ai < 1400Mpc//i or Bq < 0.2 at 
the 95% c.l. 

We have then studied theories with a Yukawa-type 
interaction between dark matter particles; these models 
require two further parameters: the coupling /3i and the 
time evolution of the scalaron mass s. In this case we 
have shown how the constraints on the theory lengthscale 
weakens, but we obtain a bound on the coupling 0.75 < 
Pi < 1.25 at the 95% c.l. when using all the data sets. 

Our result proves once more that tests of structure 
formation and attempts to reconstruct the evolution of 
the gravitational potentials are crucial to distinguish be- 
tween modifications of gravity and dark energy theories. 
Future improvements along these lines will be likely 
achieved by combining all the different available probes, 
such as the ISW, weak lensing, galaxy cluster counts 
and peculiar velocities (see for instance [H, [56l - l59| for 
forecasts) (gi} . 



VI. CONCLUSIONS 

In this paper we have presented an analysis of 
parametrised modified gravity theories in light of recent 
data from the CMB and the large scale structure of the 
Universe. 

We have focused first on the f{R) class of modified 
gravities, which are described by a single extra param- 
eter related to the mass of the scalaron, which can be 
parametrised as its present wavelength Bq . We have per- 
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expected posterior on Ai 
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1 
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baseline 


Ai < 1400 Mpc//i 


A? 


oc Ai 


> 


positive tilt 


biased to higher Ai 


Ai < 1900 Mpc//i 


Bo 


oc Ai 


> 


positive tilt 


biased to higher Ai 


Ai < 1900 Mpc//i 


Log xi 


oc 1/Ai 


< 


negative tilt 


biased to lower Ai 


Ai < 700 Mpc/h 
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deflned by transformations whose second derivative is positive (negative) will produce effective tilted priors on Ai. This will 
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